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Magnetic field-driven phase transitions have long been 
of interest in studies of magnetoresistance^ and mag- 
netic shape memoryP Recently, research into the mag- 
netocaloric effect (MCE) in the vicinity of such transi- 
tions has been revived, fuelled by interest in solid state, 
gas-free methods of cooling. 3 4 Most works that exam- 
ine materials with large MCEs focus on Curie transitions 
in ferromagnets such as Gd 5 (Si,Ge) j^, MnFe(P,Z)P and 
La(Fe,Si)i3pl They exhibit a conventional MCE, namely 
a positive change of temperature when an increasing 
magnetic field is applied. However, a smaller set of in- 
verse magnetocaloric materials have also attracted inter- 
est. In these, the MCE is associated with a field-induced 
transition to a high magnetisation state that exists above, 
rather than below, a critical temperature. As a result the 
magnetic field causes a decrease in temperature. Exam- 
ples include Hcuslers in which re-entrant ferromagnetism 
appears with increasing temperature due to the the pres- 
ence of a structural transformation to a phase with ele- 
vated Curie temperature^ and metamagnetic antiferro- 
magnets such as MnsGaCP 1 , Fe-Rh^and CoMnSP 1 , the 
subject of this article. 

Whether the MCE is conventional or inverse, its room 
temperature magnitude is normally enhanced when there 
is significant magneto-elastic coupling. This coupling 
brings about a first order magnetic transition that re- 
leases a significant fraction of the available entropy over 
a narrow temperature window. Of the inverse magne- 
tocalorics, Feo.49Rho.51 still holds the record for the mag- 
nitude of the adiabatic temperature change per Tesla 
of applied fielcPH due to its close-to-optimal value of 
dTt/dH, the rate at which the metamagnetic transi- 
tion temperature changes with applied fieldPSi An in- 
vestigation of magneto-elastic coupling in inverse mag- 
netocaloric materials is therefore expected to shed light 
on mechanisms for achieving large MCEs, partly through 
tuning dTt/dH to more optimal values than are seen in 
most ferromagnetic materials. 

Recently we used high resolution neutron scattering to 



observe giant magneto-elastic coupling in CoMnSi and 
Coo.gsNio.osMnSi, associated with the temperature evo- 
lution of the antiferromagnetic state of both materials. 
In each, the nearest neighbour Mn-Mn separations were 
shown to change by about 2% over a 150 K range, the 
largest such change seen in a metallic magnet.^ Further- 
more, application of a magnetic field suppresses the heli- 
cal antiferromagnetism and yields a metamagnetic tran- 
sition to a high magnetisation state. For temperatures 
close to the Neel temperature, this transition is continu- 
ous. At lower temperatures the metamagnetic transition 
couples with the underlying change in Mn-Mn separa- 
tion and goes through a tricritical point to become first 
order, accompanied by an enhanced inverse MCE. Such 
tricriticality is a useful property, combining large entropy 
changes with low hysteresis. Indeed all of the first order 
magnetic refrigerants (MnFeP-based, La(Fe,Si)i3-based 
and manganites) currently trialled in prototype magnetic 
cooling engines exhibit field-induced critical points. 

Since metamagnetism is known to appear with dop- 
ing as an intermediate state between ferromagnetism and 
antiferromagnetism in other Prima Mn-based alloys, our 
findings in CoMnSi led us to a simple theoretical model 
for the existence of either ferromagnetic or antiferromag- 
netic groundstates in such Mn-based structures, based 
on how close their nearest neighbour Mn-Mn separations 
are to a critical value of around 3 AP^ That theoret- 
ical work has also been widened to the prediction and 
synthesis of new CoMn(P,Ge) metamagnets from ferro- 
magnetic end-members, by tuning the nearest neighbour 
Mn-Mn distance towards the same critical valued 

In this article we extend our analysis of giant magneto- 
elastic coupling in CoMnSi to a series of substituted ma- 
terials in which local exchange interactions, and hence 
metamagnetism, are altered. We demonstrate the generic 
presence of giant magneto-elastic coupling in all com- 
pounds studied, and how its magnitude affects the sensi- 
tivity of the antiferromagnetic state to an applied mag- 
netic field, and hence tricriticality. Furthermore, we ex- 
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amine the nature of the metamagnetic transition and 
its entropic constitution by analysing heat capacity and 
magnetostriction data. We show that there is a large 
change of electronic entropy at the transition that is 
partially compensated by a smaller change in lattice en- 
tropy and that the entropic balance corroborates our den- 
sity functional model of the V-shaped bandstructure of 
CoMnSi alloys near to the Fermi energy. 

The remainder of this article is organised as follows: 
section |T] gives details of our synthesis and characterisa- 
tion work. Results are given in Section |n] and discussed 
in Section Hill 



I. EXPERIMENTAL 

The samples that we study here were formed by co- 
melting appropriate amounts of high purity elements Co 
(99.95 %), Mn (99.99 %), Ni (99.994 %), Fe (99.995 %), 
Cr (99.995 %) and Si (99.9999 %) under argon in an 
induction furnace, followed by post-annealing and slow 
cooling. Details are given elsewhereP^l An oxide layer 
on the as-received manganese was removed by chemical 
etching.^ In order to avoid potential embrittlement by a 
martensitic transition between the orthorhombic ground- 
state and a high temperature hexagonal state, the alloy 
was cooled slowly during solidification. The temperature 
of the martensitic transition was established by simul- 
taneous differential thermal analysis (SDT) in a TA In- 
struments Q600. After each alloy ingot was first formed 
it was wrapped in tantalum foil, vacuum sealed in a sil- 
ica tube, and annealed for 60 hours at either 1123 K 
or 1223 K to fully relieve strain, before cooling to room 
temperature at 0.2 Kmin -1 . 

Magnetic measurements in fields of up to 9 Tesla 
were performed on poly crystalline samples with a Cryo- 
genic Ltd. vibrating sample magnetometer (VSM) and 
a Quantum Design PPMS VSM. AC Susceptibility mea- 
surements in a Lake Shore 7130 susceptometer were per- 
formed over a range of frequencies and between 20 and 
320 K. The Neel transition temperature of several alloys 
was measured in zero magnetic field in a conventional 
heat flux differential scanning calorimeter (DSC, TA In- 
struments Q2000) with a temperature range of 120 K to 
873 K. Small amounts of alloy (~5 mg) were encapsulated 
in an aluminium pan and an empty aluminium pan was 
used as a reference sample. The heating rate in the DSC 
was set to 10 Kmin -1 and data was recorded between 
273 K and 673 K. Heat capacity measurements were also 
employed in a separate apparatus to determine the Debye 
temperature and the Sommcrfcld coefficient. That appa- 
ratus was designed for experiments at low temperatures 
under quasi-adiabatic conditions. A detailed descrip- 
tion of the experimental setup is given elsewhere.^ The 
calorimeter consisted of an adiabatically shielded sample 
holder equipped with a resistive heater (a strain gauge) 
and a thermometer (a Cernox sensor). Quasi-adiabatic 
conditions, and hence accurate operation of the calorime- 



ter were limited to temperatures between 1.5 K and 180 
K. 

Structural characterisation was carried out by X-ray 
and neutron diffraction. We conducted a Rietveld refine- 
ment of data from room temperature X-ray diffraction 
using Cu-Ka radiation. X-ray diffraction at room tem- 
perature showed that most samples were single phase 
to within experimental resolution. Where not, this is 
indicated in Table [l] As in our previous worlP^, struc- 
tural refinements of these data agreed well with neutron 
diffraction data results close to room temperature. We 
therefore only present unit cell refinements obtained by 
high resolution neutron diffraction. 

Neutron diffraction was carried out at the time-of- 
flight High Resolution Powder Diffractometer (HRPD) 
and at the General Materials Diffractometer (GEM) at 
ISIS, UK. The former has a resolution of Ad/d - 1 x 10~ 4 
and was used at temperatures between 4.2 K and 500 K. 
The latter is capable of hosting a 7 Tesla magnetic field 
and was used for magnetostriction studies of a CoMnSi 
ingot between 150 K and 300 K. The choice of a single 
ingot was to limit sample movement in the field. The 
ingot was held at the bottom of a vanadium can by a 
rolled sheet of cadmium. Magnetic field steps were cho- 
sen based on magnetisation measurements. The incom- 
mensurate magnetic structure and its variation with tem- 
perature and applied magnetic field are difficult to resolve 
accurately without a single crystal sample and will not 
be discussed in detail here. 



II. RESULTS 
A. Structure and magnetism 

The Neel transition in zero magnetic field is contin- 
uous, as evidenced by a broad peak in heat capacity 
measurements that occurs at T/v~380 K for all sam- 
ples (Figure [I]). Despite the small difference in Neel 
temperature between the alloys, there is a strong vari- 
ation in thermal expansion behaviour with composi- 
tion. Figure [2] shows the zero-field thermal expan- 
sion of CoMnSi, Coo.gsNio.osMnSi, CoMno.9sCro.02Si and 
CoMno.95Feo.05Si from HRPD data. These are the four 
main samples involved in this study. As in our previous 
studjEl, diffraction data were first recorded at the lowest 
temperature, before heating the sample to the other tem- 
peratures at which data was taken. If additional patterns 
were collected the sample was first cooled down to base 
temperature and then heated up to the desired tempera- 
ture in order to eliminate any possible effects of thermal 
hysteresis. At each step the temperature was equilibrated 
for 20 min and a total neutron current of either 75 fiA or 
60 jLtA was collected for each frame. We used the GSAS 
code^ to perform a Rietveld refinement of the crystal 
structure. Room temperature structural parameters are 
shown in Table |U 

Figure [3] shows the temperature evolution of the two 
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FIG. 1: (Color online) Specific heat capacity as a function 
of temperature for four CoMnSi-based alloys. The anomaly 
around 380 K is assigned to the Neel phase transition, Tm- 
Measurements were conducted in zero magnetic field. De- 
spite the small difference in Neel temperature (shown by the 
arrows), the alloys have quite different magnetic characteris- 
tics and strong variations in thermal expansion. Symbols and 
colours are the same as in Figs. [2] and [3] 
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FIG. 2: (Color online) Temperature dependence of lattice pa- 
rameters measured on powders using HPRD. Similar features 
are seen in all 4 samples. The large change in cell parameters, 
including an a axis NTE, occurs well below the Neel transi- 
tion temperature, TV. Such features are strongest in the Ni- 
and Fe-substituted samples. 



nearest neighbor Mn-Mn distances, di and d^ in zero 
field, obtained from Rietveld refinement of the same data. 
d\ and d-^ are the distances between different chains of 
manganese atoms, and between Mn atoms on the same 
chain, respectively.^ We previously identified d\ as the 
key structural parameter that determines whether Mn- 
based orthorhombic magnets in the Pnma space group 
have a f errom agnetic ground state or an antiferromag- 
netic on d 12 * 15 *, since, in a certain range of di the coupling 
between Mn-Mn atoms on different chains is antiferro- 
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FIG. 3: (Color online) Manganese nearest-neighbor distances 
as a function of temperature for the same four samples as in 
Fig. [2] Distance di is between different chains of manganese 
and d2 is between Mn on the same chain. 



magnetic. The temperature variation of d\ in differently 
substituted compounds is therefore of interest here. 

The structural data for CoMno.9sCro.02Si and 
CoMno.95Feo.05Si allow us to compare samples in which 
the helical antifcrromagnetism is strengthened or su- 
pressed, as can be concluded from the magnetisation 
data, which we now present. Figure [4] compares isother- 
mal magnetisation curves at 300, 250 and 200 K in fields 
of up to 8 Tesla for the same four samples as in Figure [2] 
We see that, relative to the stoichiometric material, Cr- 
substitution increases the critical fields. Tricritical points 
are observed in all the alloys, and can be seen from the 
onset of hysteresis in the magnetisation curves in Fig. [4] 
and the critical field vs temperature phase diagram that 
can be constructed from them. Figure [5] shows the tem- 
perature variation of magnetic susceptibility measured at 
137 Hz in a 100 A/m field. The Fe- and Co-substituted 
alloys have approximately 10 times greater initial suscep- 
tibility than CoMnSi, indicating weakened antifcrromag- 
netism (AFM) and greater proximity to ferromagnctism 
(FM). 

Considering Figs.[T][2j and [4] together we see that the 
features of negative expansion along a, and positive ex- 
pansion along b and c are found in all alloys, as well as the 
giant magneto-elasticity of between 1% and 2% change 
in Mn-Mn separation with temperature, well below TV. 
However, the largest changes in a and in d\ are seen in 
the Ni- and Fe-substituted materials, which also have the 
lowest critical fields. Cr would seem to increase the crit- 
ical fields and suppress the magneto-elastic interactions 
that lead to metamagnetism and tricriticality. 



B. Entropy change and magnetostriction 

Our interest in tricriticality stems from a desire to un- 
derstand, control and tune the feedback between struc- 
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sample ID 


cms38-a 


cms40-a 


cnms39-a 


cmfs33-a 


cmcs41-a 


formula 


CoMnSi 


CoMnSi 


Co .95Ni .o5MnSi CoMno.95Feo.05Si 


CoMno.98Cro.02Si 


T anncal (K) / t (h) 1223 / 60 1223 / 60 


1223 / 60 


1123 / 60 


1223 / 60 


2 nd phase 


t 


u 


11 


y 


11 


a (A) 


5.8666(1) 


5.8677(2) 


5.8527(2) 


5.8398(2) 


5.8725(5) 


b(A) 


3.6862(9) 


3.6865(1) 


3.6852(1) 


3.6924(2) 


3.6820(3) 


c(A) 


6.8517(2) 


6.8528(3) 


6.8618(2) 


6.8607(3) 


6.8468(6) 


v (A 3 ) 


148.1724 


148.2352 


147.9967 


147.9382 


148.0453 


XCo 


0.1607(6) 


0.1544 


0.1550(5) 


0.1614(7) 


0.152(1) 


ZCo 


0.0605(5) 


0.0483 


0.0540(5) 


0.5563(6) 


0.5563(9) 


XMn 


0.0146(7) 


0.02225 


0.0204(7) 


0.0244(8) 


0.023(1) 


ZMn 


0.6816(4) 


0.6774 


0.6764(4) 


0.1884(6) 


0.1793(8) 


XSi 


0.269(1) 


0.2664 


0.261(1) 


0.264(2) 


0.767(2) 


ZSi 


0.3825(9) 


0.3595 


0.3744(8) 


0.378(1) 


0.617(1) 


flaw 


20.8 


20.5 


18.0 


17.8 


26.0 


^struct 


119(P 


119(P 
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TABLE I: Table of room temperature structural data, derived from Rietveld refinement of HRPD neutron diffraction data. 
Also shown are: annealing conditions (Tanneai and duration, t); presence of a 2 nd (hexagonal) minority phase (yes/no/trace); 
goodness of the Rietveld refinement (Rnw) and the high temperature martensitic phase transition temperature T str uct (where 
known) . 



ture and magnetism which seem to be key to the gi- 
ant magneto-elasticity and large magnetocaloric effect 
found in these materials. We focus on the temperature 
range in which the magnetic field induces a metamagnetic 
transition, below T/v, with an associated inverse magne- 
tocaloric effect. (There is a smaller conventional MCE 
associated with the depolarisation above T/v). Isother- 
mal magnetisation data presented above are used to call- 
culate isothermal entropy change in each material, in a 
given field change. We use the Maxwell relation: 

ASiT^^f^^M, (1) 

which can subtlely overestimate the equilibrium phase 
change entropy around hysteretic phase transitions.^ 
Such signs of irreversibility are measurable in CoMnSP^l 
but are not the subject of the current study. Instead 
we use the Maxwell-derived entropy change to compare 
materials since thermomagnetic hysteresis is only evident 
above the tricritical field. 

Figure [6] shows the derived isothermal entropy changes 
for all four alloys, in fields of up to 5 Tesla (for CoMnSi 
and Corj.95Nio.o5MnSi) or 9 Tesla for the two alloys most 
magnetically distinct from CoMnSi: CoMno.95Feo.05Si 
and CoMno.98Cro.02Si. As can be expected from the 
M (H) plots, there is a low-field enhancement of the en- 
tropy change seen in the Fe-and Ni-substituted materials 
when compared with CoMnSi. The Fe-substituted case 
was presented previously^, with a focus on the entropy 
change measured directly by calorimctric methods on a 
100 /im-sized fragment. The additional data taken here 
enables us to clarify a statement made there. On increas- 
ing the maximum magnetic field strength the metamag- 



netic phase transition temperature decreases monotoni- 
cally (see Figure [4]) while the entropy change found at 
the transition first appears to increase, reaches a maxi- 
mum, and then becomes smaller. However, the appar- 
ent downturn in the entropy change is an analysis arte- 
fact when a hysteretic phase transition is broadened, and 
hence incomplete, due to effects such as polycrystalline 
inter-grain strain. As an example, in a field change of 
9 Tesla the lowest temperature at which the phase transi- 
tion is complete in an increasing field in CoMno.95Feo.05Si 
is 95 K, which in fact corresponds to the transition tem- 
perature for a 8.5 Tesla field. In addition the magneti- 
zation jump and the gradient of the magnetic phase line 
in (iJ, T) space are both approximately constant from 
3 Tesla < ^.qH < 9 Tesla, so from a Clausius-Clapeyron 
analysis a significant downturn in entropy change is not 
expected. 

Turning to the contributions to the entropy change, we 
previously concluded that any change in lattice entropy is 
negative and so adds in opposition to any change in mag- 
netic or electronic entropypS However, the signal used to 
provide a measure of volume change in those experiments 
was the change in resistance of the calorimeter mem- 
brane. Since the thermal expansion of CoMnSi-based 
materials is highly aniostropic and the fragment was very 
small we could not be certain of the size or sign of the 
volume change in our materials. Likewise, in previous 
capacitance dilatometry studies of Coo.gsNio.osMnSp^, 
textured polycrystals were used and so the true depen- 
dence of sample volume on magnetic field was unknown. 
Therefore, we here examine the magnetostriction of our 
materials directly using neutron diffraction and establish 
the origin of large, competing, changes in entropy that 
constitute the metamagnetic phase transition. 
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FIG. 4: (Color online) Upper plot: critical field vs tempera- 
ture phase diagram for the metamagnetic transition of the 4 
main alloys studied here. Lower 3 plots: isothermal magneti- 
sation data at 300 K, 250 K and 200 K as examples of the 
data used to construct the phase diagram. Tricritical points 
are seen from the narrowing of the hysteresis. Data for the 
cms40-a sample of CoMnSi are shown. 



Figure [7] shows the transverse magnetostriction in a, b, 
c, and volume for CoMnSi ameasured on GEM. CoMnSi 
diffraction patterns were collected at 300, 285, 270, 255 
and 225 K during cooling with a total counting time of 
around 30 minutes per pattern. We see from Fig. [7] that 
magnetostriction has similar properties to thermal ex- 
pansion: namely that a-axis magnetostriction is negative 
while that along b and c is positive. The resulting volume 
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FIG. 5: (Color online) Magnetic susceptibility of the 4 main 
alloys studied here measured at 137 Hz in a 100 A/m field. 
The susceptibility of CoMnSi and CoMno.9sCro.02Si is lower 
than the other two alloys, and are plotted on a separate axis. 
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FIG. 6: (Color online) Isothermal entropy change of the 4 
main alloys studied in integer field steps up to 5 Tesla (up- 
per plots) or 9 Tesla (lower plots). We see that increased 
proximity to a room temperature tricritical point in Fe- and 
Ni-doped samples increases the low field entropy change. The 
color scheme used is the same in all 4 sub-plots, and up- or 
down-triangles represent the effect of applying the Maxwell 
relation to increasing-field or decreasing-field M(H) data. 



change is small in fields of up to 3 Tesla, above which the 
a-axis behaviour dominates and a large negative mag- 
netostriction is observed. In the highest applied mag- 
netic fields the lattice parameters assume values never 
reached during thermal expansion measurements of Fig- 
ure^ For example, in a magnetic field of 6 Tesla the ex- 
axis shrinks almost twice as much as it does on changing 
the temperature from 4 K to 330 K. On the other hand, 
the deformation along the 6-and c-axes induced by the 
same magnetic field change are comparable to the effect 
of temperature. This suggests that the magnetostriction 
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FIG. 7: Magnetic field dependence of lattice parameters of 
CoMnSi (cms40-a) as a function of temperature. The zero 
field values agree well with the data obtained from HRPD on 
a different sample of CoMnSi (cms38-a). The magnitude of 
the lattice parameter changes increases with decreasing tem- 
perature. 
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FIG. 8: Low temperature specific heat data for four CoMnSi- 
based metamagnets: the C/T axis intercept is the Sommerfeld 
coefficient, 7, and is tabulated in Table |TT| 



is influenced predominantly by changes in magnetic order 
along the a- axis. 

The field-induced change in lattice entropy, as dis- 
cussed in section [Hi] below, is thus expected to be large 
and negative above ^3 Tesla. The origin of the positive 
entropy change that dominates the overall MCE is now 
presented by examining the electronic degrees of free- 
dom via the Sommerfeld coefficient, 7. In the meta- 
magnet Fe-Rh two strategies were previously adopted 
to examine 7, both using low temperature heat capac- 
ity measurements. The first approach compared the 7 
values of compositions with ferromagnetic (FM) ground- 
states to those which were antiferromagnets (AFM)P^ 
The FMs were found to have significantly higher 7 val- 
ues: A7 = Jfm ~ Iafm ^6-7 mJ/mol K 2 . A similar 
value of A7 ~ 14 mJ/mol K 2 was found for a lower- 
temperature metamagnet, CeFeo.9Coo.1PSl The second 
approach examined the magnetic field-induced change in 
7 in metamagnetic compositions of (Fei_ a Ni x )o.49BJio.5i 
with x ~ 0.03™ Again a similar increase in 7 was seen 
on entering the high magnetisation state, thus leading to 
the conclusion that the giant MCE at the metamagnetic 
transition was due in large part to an increase in the elec- 
tronic degrees of freedom, and therefore to a change in 
the density of states at the Fermi energy. This increased 
density of states was a hallmark of the high magnetisa- 
tion state, brought about either by chemical substitution 
or by a magnetic field. 

We employ the first approach here since the fields 
required to induce metamagnetism at very low tem- 
peratures in CoMnSi are much higher than those in 
Fe-Rh. In Figure [8] we show specific heat data on 
four different CoMnSi-based compositions. Antifcrro- 
magnetic compositions are represented by CoMnSi and 
CoMnSio.95Geo.05; the synthesis and magnetic charac- 
terisation of the latter composition was described pre- 



viously.^ In order to compare a CoMnSi-based material 
with a ferromagnetic groundstate, a composition is re- 
quired that has a higher level of chemical substitution 
than is present in any of the samples presented in the 
above magnetisation and neutron experiments. There- 
fore two new compositions were prepared for this experi- 
ment using the same protocol as given in Section]!] These 
are Coo.5Nio.5Mno.sCro.2Si and Coo.5Nio.5Mno.9Cro.1Si. 
They were found to have high magnetisation ground- 
states (not shown here) which are assumed to be fer- 
romagnetic. 

We can see from Figure [8] that, as with Fe-Rh, there is 
a very large increase in 7 when comparing a FM compo- 
sition with a AFM one. Compared to Fe-Rh, the change 
in 7 in CoMnSi-based materials would seem to be much 
larger, at around 15 mJ/mol K 2 . 



III. DISCUSSION 

A. Giant magneto-elasticity and tricriticality 

The results presented here enable us to establish and 
quantify several properties of CoMnSi-based metamag- 
nets for the first time. From section III Al we see that 
the effects of Fe or Ni doping are opposite to that of Cr. 
Critical metamagnetic fields are reduced be either of the 
former, and enhanced be the latter. The peak in heat 
capacity at Tjv is suppressed by Fe or Ni, and enhanced 
by Cr. The change in lattice parameters and in Mn-Mn 
distances (in particular d%) become stronger with Fe, are 
shifted down in temperature by Ni, but are suppressed by 
Cr. If we ascribe the strength of heat capacity jump, size 
of metamagnetic critical field and decreased magneto- 
elastic coupling to the presence of antiferromagnetism 
with increased anisotropy then, as we now outline, the 
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above data can be viewed as self-consistent. 

In CoMnSi, the incommensurate helical magnetic or- 
dering wavevector is (0 q), where q ~ 0.4 at low tem- 
perature. One possibility is that Fe or Ni doping reduce 
the value of q, bringing the material closer to ferromag- 
netism and reducing the metamagnetic critical field. Cr 
would seem to have the opposite effect of stiffening the 
helical arrangement of moments, presumably as a result 
of strengthening the antiferromagnetic exchange that is 
mediated by the inter-chain Mn-Mn separation d\. We 
previously proposed a simple model of magnetic ground- 
states in Mn-based Prima alloys based on MnP, where low 
values of d% suppress any magnetic order, and then, on 
increasing d\ , two boundaries between different magnetic 
groundstates are seen. First FM, then AFM and finally 
FM states are stabilised as d\ increases.^ The d\ val- 
ues at the theoretical FM/AFM phase boundaries^ are 
~3 A and ^3.36 A. The compounds presented in Figs. [2] 
and [3] have d\ values between 3.06 and 3.08 A at 4 K 
and are therefore clearly on the lower boundary betwen 
low volume ferromagnetism and higher volume antifcr- 
romagnetism. CoMno.9sCro.02Si has the highest value 
of d\ at low temperature and can therefore be thought 
of as sitting furthest from that magnetic phase bound- 
ary. The findings of high critical fields and reduced 
magneto-elastic coupling in this material are therefore 
self-consistent. 

These conclusions are also in agreement with the ob- 
servation that Cr reduces the magneto-elastic coupling in 
zero field, if indeed that coupling arises from a difference 
between the value of d\ in the antiferromagnetic ground- 
state and in the high magnetisation state and that the 
reverse is true in the case of the Fe and Ni substitutions 
studied. Another interesting possibility is that the Mn 
moments, which are believed to lie in the (a,b) plane in 
CoMnSi, are tilted out of plane to form a canted fer- 
romagnetic structure, thus enhancing the magnetic sus- 
ceptibility and reducing the (tri)critical fields. Further 
neutron diffraction measurements and single crystalline 
samples are clearly required to examine the relationship 
between magnetic and crystal structure. Another study 
of the magnetoresistance of CoMnSi has suggested the 
presence of a spin reorientation transition below the Neel 
temperature on the basis of a cusp in the ac susceptibility 
at -150 K. 27 However, our magnetic susceptibility mea- 
surements (Figure [5| show no sign of such a feature at 
this temperature. 



B. Electronic contributions to AS 

We now present an analysis of the components of the 
isothermal entropy change, by using as an example the 
transition in CoMnSi at about 225 K and in ~5 Tesla, 
where the change of volume is around —0.3% (Fig- 
ure [7|) AS is about 7 J kg" 1 K _1 . In order to deconstruct 
the entropy change at this transition we start with the 
electronic contribution and then proceed to the role of 



composition 7 (mjmol *K 2 ) 7 (mjkg X K 2 ) 



CoMnSi 


5.2 


36.6 


CoMnSio.95Geo.05 


6.7 


46.5 


Coo.5Nio.5Mno.9Cro.1Si 


19.2 


135.6 


Coo.5Nio.5Mno.8Cro.2Si 


21.3 


150.8 



TABLE II: Values of 7-coefHcient determined from a linear fit 
of the low temperature specific heat capacity data. 



the lattice. It is not straightforward to calculate mag- 
netic entropy contributions without detailed knowledge 
about the magnetic structure in both magnetic states and 
whether there is a jump in the site moment as found in 
other systems such as FeRh and MnFe(P,Z)PHSI With- 
out any evidence for such changes in site moment, mag- 
netic contributions to the phase transition between two 
ordered states are expected to be small and will not be 
considered other than through a change in the density 
of electronic states (DOS). Hereafter we refer simply to 
the "electronic" entropy change, while recognising that 
DOS changes will contribute to the itinerant part of the 
magnetic entropy change. 

The large change of electronic entropy at the metamag- 
netic transition is consistent with our earlier theoretical 
result that a helical magnetic state has a lower density of 
electronic states at the Fermi energy than a (fictitious) 
ferromagnetic ground stated It is also perhaps consis- 
tent with the extraordinary sensitivity of CoMnSi to ex- 
ternal parameters such as synthesis conditions and ex- 
ternal parameters such as pressure and magnetic field. 
We previously showed a variation of 200 K in the low 
field metamagnetic transition temperature from litera- 
ture dataP^ Even in the work presented here, two ingots 
of the same nominal composition (cms38-a and cms40-a) 
have slightly different critical fields and temperatures. 

A list of the 7- coefficients extracted from the data in 
Figure [T] is given in Table |H| It shows that the electronic 
specific heat capacity which is linearly proportional to 7 
is approximately tripled when moving from an antifer- 
romagnetic alloy to a ferromagnetic one. It can reason- 
ably be assumed that the electronic heat capacity at the 
field-induced metamagnetic phase transition in CoMnSi 
changes by a similar order of magnitude. 

Although this assumption is reasonable, it is oversim- 
plified for two reasons, the impact of which need to be 
examined briefly. Firstly, the substitution used to gen- 
erate ferromagnetism is not isoelectronic. This makes a 
direct comparison of the 7-coefficients of two compounds 
nontrivial. Secondly, even if such an experiment can 
be conducted with isoelectronic samples it can still not 
be excluded that some change of 7 comes from volume 
changes and their impact on the density of states at the 
Fermi energy. However, this effect scales as the fractional 
volume changed AS^volume = 1^7^ an d is therefore 
small (~ 0.016 J kg" 1 K" 1 at 225 K). Therefore the three- 
fold increase in 7 observed is unlikely to originate from 
substitution only, or from volume changes. 



A band structure feature is therefore most likely to 
contribute to such a large difference between the entropy 
of the AFM state and field-induced high magnetisation 
state. This idea is corroborated by the large difference 
between the DOS of a non-collinear AFM and that of 
a collinear FM structure in CoMnSi that we previously 
calculated from DFTP^ However, the shape of the DOS 
is crucial here. It might seem simplest to examine the 
average difference between the electronic heat capacity 
in the AFM and FM alloys which is -400 mJkg" 1 K~ 2 
and then make an estimate of the electro nic en tropy 
change, as is often done in other material d 25 * 26 ! using 
A5 c i = TtA'y. In the current case, this yields a value 
of 23 Jkg~ 1 K _1 at 225 K which would be an upper 
limit for AS e \ at this temperature. However, the form 
of AS C \ oc T t predicts that the electronic contribution 
should grow linearly in magnitude with T t (decreasing 
field) which contradicts the opposite trend in AS(T t ) 
seen in Figure [6] for most of the field range studied. In 
other words, the effects of finite temperature on energy- 
dependent bandstructure need to be taken into account. 

We can see from our earlier DFT calculations® (Fig- 
ure 4 of that study) that the AFM DOS has a dip near 
the Fermi Energy, Ep which has an energy width of order 
ksT and that the lowest DOS value is only to be found 
close to Ep itself. Thermal occupation of the DOS is 
indeed therefore relevant to the measured electronic en- 
tropy change. We may hence use a simple model of a 
V-shaped DOS, of width 0.1 eV and employ the statis- 
tical description of fcrmionic entropy to the FM and the 
AFM DOS: 

Sd = -k B J d£[/log/ + (l-/)log(l-/)]. (2) 

The (dis)occupation of states with DOS higher than the 
dip at the Fermi energy reduces the entropy difference 
between FM and AFM states, relative to the simplified 
model stated above by a factor of two or more in the room 
temperature range (see Figure l9j). The precise reduction 
depends on the gradient of the DOS with energy on either 
side of Ep. Parameters here are chosen to match the form 
of the DOS calculated previously and the approximate 
increase in DOS at the Fermi level obtained from that 
calculation and from Tabic Ull 

Hence, rather than 23 Jkg _1 K _1 at 225 K, a value 
closer to 10 J kg -1 K _1 electronic entropy change would 
seem reasonable at this temperature. A fuller examina- 
tion of the entropy change across a range of tempera- 
tures — including at low fields where the transition to 
a high moment state is incomplete — would perhaps re- 
quire more than the V-shaped AFM-DOS model. The 
motivation here is to stress the importance of the DOS 
profile near the Fermi Energy. Despite the appearance of 
a very large change in 7 coefficient, the accessible change 
in electronic entropy is limited by the varation of AFM 
DOS near to Ep . 

In order to confirm the large electronic entropy change 
in CoMnSi further experiments on isoelectronically sub- 




Temperature (K) 

FIG. 9: A simple V-shaped model DOS for the AFM state 
near Ef (red solid line, inset) yields a reduced entropy change 
(red solid line, main figure) relative to that of two energy- 
independent DOS in both the AFM and FM states) . 

stituted compounds are encouraged. A possible alloy se- 
ries could be germanium rich CoMnSii-zGez which are 
ferromagnetic for large x!^ By changing the germanium 
content it may be possible to tailor the magnetic prop- 
erties of the alloy in order to access the transition in the 
magnetic field of a laboratory magnet. 

C. Lattice contributions to AS 

The lattice contributions to AS can be estimated us- 
ing the protocol described by Jia et alPSl Accordingly, 
there are two components to the lattice entropy: (i) a 
phonon contribution due to the shift in Debye tempera- 
ture caused by the field and (ii) an elastic contribution 
due to the deformation of the bulk material. The De- 
bye temperature, Td can be estimated from a fit of the 
experimental heat capacity data to the a Debye model 
for GV(T). This assumes that Cp(T) data taken under 
isobaric conditions differ only slightly from the isochoric 
data required for the Debye model. The difference be- 
tween isochoric and isobaric C(T) requires knowledge of 
the isothermal compressibility kt and the coefficient of 
linear thermal expansion, a^, of the material: 

C P -CV = ^ (3) 

where v is the molar volume. Due to the small thermal 
expansion of CoMnSi at low temperatures, this correction 
does not contribute significantly to the quality of the fit, 
especially since the model is a rather simple one. 

A Debye temperature in the antiferromagnetic state of 
Td = 410 K results in the best fit of experimental data to 
the Debye model for temperatures T < 200 K. Strictly 
speaking the Debye model is only valid for monatomic 
isotropic materials. In CoMnSi the lattice expands 
anisotropically which might explain the deviation from 
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the simple model at temperatures above this range. In- 
terestingly the temperature range above which the agree- 
ment of experimental heat capacity data with the model 
starts to get worse coincides with the start of negative 
thermal expansion along the a-axis (Figure [2]) . 

The lattice entropy for a system is then given by 33 34 : 



SiAT) = -3NRln(l - e 



+12NR 



T 



T d /T 



1 



-dx 



(4) 
(5) 



where N is the number of atoms per mole and R the gas 
constant. A volume change influences the phonon system 
and hence the Debye temperature by23 



AV 



T D 77 V 



(6) 



where rj is the Griineisen parameter. This parameter 
depends on the material and is between unity and 3 for 
many systems. In LaFei3^ x Si x systems r\ has elsewhere 
been set to 6 which resulted in the best fit of experimental 
data to a model similar to the one discussed herePS The 
Griineisen parameter for CoMnSi is unknown. Setting 
j) = 3 as an example yields a change in Td of about 3.7 K 
or a x^ lshM = 413.7 K in the spin-aligned state. The 
lattice entropy change can then be calculated from the 
difference between 5j a t values at a single temperature 
(225 K) using two different Debye temperatures. For 
i] = 3 we obtain A5i a t — — 4Jkg _1 K _1 . The sign of the 
lattice entropy change is negative since a higher Debye 
temperature in the spin-aligned state means that fewer 
phonon modes are excited at a fixed temperature. 

The elastic contribution, A5 e i a examines the deforma- 
tion of the crystal as a whole. According to Hooke's law 
the energy difference between two states assuming that 
the bulk modulus is constant is given by: 



where Vb is the reference volume. For CoMnSi at 225 K 
the volume difference between the AFM and FM-like 
states is AV = 7.6 x 10 -2 cm 3 /mol. The unit cell volume 
in the AFM reference state at 225 K is 22.342 cm 3 /mol. 
The bulk modulus of CoMnSi is unknown but we might 
assume that it is equal to that of MnSi, B = 1.63 x 
10 7 N/cm 2 !^ From these values we can obtain a tiny ef- 
fective latent heat at 225 K of A,S* c i a = 0.66 Jkg" 1 K -1 . 
It should be emphasised that this value is based on an es- 
timated, field- and temperature-independent bulk mod- 
ulus. 

In summary it seems that in CoMnSi the positive field- 
induced isothermal entropy change is mainly due to an 
contribution from the electronic density of states (similar 
to FeRh-based alloys which also exhibit an antiferro- to 



ferromagnetic phase transition and estimated here to be 
10 Jkg _1 K _1 at 225 K) and counteracting phonon den- 
sity contributions (e.g. around -4 Jkg _1 K _1 at 225 K). A 
naive isotropic analysis of the phonon entropy contribu- 
tion is clearly not sufficient. Direct phonon spectroscopy 
is needed — as it is for several magnetocaloric systems — 
in order to more accurately determine the role of phonons 
in contributing to the total entropy change determined 
from magnetic or calorimetric measurements. The role of 
the lattice can be extremely important in other materials, 
either in adding to the magnetic entropy change (Gd-Si- 
Ge^3 or in counteracting it, as seen here and in La-Fe- 
SiP2l Recent theoretical work by Basso has extended the 
model of first order phase transitions by Bean and Rod- 
bell to calculate the role of the lattice in the isothermal 
entropy change at a magnetic field-driven phase transi- 
tion! 3 -^ Further extension of such a model to anisotropic 
systems such as found in MnFe(P,Z) and CoMnSi will 
be of interest in order to relate anisotropy in thermal 
expansion to magneto-elastic entropy. 



IV. CONCLUSIONS 

We have found that different elements have contrast- 
ing effects on the position of the tricritical point in 
(T, H) space and that this can be related to the size of 
the magneto-elastic coupling as manifested by the tem- 
perature dependence of manganese nearest-neighbor dis- 
tances. The entropy change associated with metamag- 
netism in CoMnSi is dominated by a very significant 
change in electronic entropy, compensated by an oppo- 
site phonon contribution. Our work motivates further 
attempts at single crystal synthesis in order to undertake 
neutron diffraction experiments to examine the temper- 
ature and field evolution of the magnetic structure, as 
well as neutron spectroscopy to examine the nature of 
the changing phonon density of states. The variation in 
tricritical fields are also an interesting test bed for finite 



g 

(7) temperature theories of magnetism! 39 * 40 ! 
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